Effects of short-range interactions on transport through quantum point contacts: 

A numerical approach 
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We study electronic transport through a quantum point contact, where the interaction between 
the electrons is approximated by a contact potential. Our numerical approach is based on the 
non-equilibrium Green function technique which is evaluated at Hartree-Fock level. We show that 
this approach allows us to reproduce relevant features of the so-called "0.7 anomaly" observed 
in the conductance at low temperatures, including the characteristic features in recent shot noise 
measurements. This is consistent with a spin-splitting interpretation of the process, and indicates 
that the "0.7 anomaly" should also be observable in transport experiments with ultracold fermionic 
atoms. 
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I. INTRODUCTION 

One of the most prominent quantum phenomena in 
mesoscopic physics is the effect of conductance quanti- 
zation. The conductance of a quantum point contact 
measured as a function of an applied gate voltage ex- 
hibits plateaus at integer multiples of the conductance 
quantum, Go = 2e'^/h, where — e is the electron charge 
and h is Planck's constant [H, i, Q. These steps are 
well understood in terms of non- interacting electrons [3| . 
But experimental conductance curves frequently show 
an additional plateau-like feature below the first con- 
ductance step at a value around 0.7 x Gq. This so- 
called "0.7 anomaly" was first investigated experimen- 
tally by Thomas et al. [H] who particularly looked at the 
magnetic field- and temperature dependence of the addi- 
tional plateau. They found that the 0.7-feature develops 
smoothly into the Zeeman spin-split plateau at 0.5 x Go 
by applying a parallel in-plane magnetic field. That is 
why those authors related this anomaly to the spin degree 
of freedom of the electrons. They conjectured the pres- 
ence of spin polarization in quasi one-dimensional junc- 
tions. In addition, Thomas and coworkers revealed that 
the 0.7 plateau becomes more pronounced if the temper- 
ature is increased. 

Since those first measurements there has been much 
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its', TI^ and theoreti- 
,[11, 24, 25, 26, 27] effort 
to explain the origin of this effect. However, a complete 
understanding is still missing. Experiments show a zero- 
bias peak in the differential conductance typical for the 
Kondo effect [loj . Furthermore, the temperature depen- 
dence can be characterized by a single parameter which 
was interpreted as the Kondo temperature. In a recent 
experiment [l^ a static spin polarization was measured, 
which contradicts the Kondo interpretation. Shot noise 
measurements [H, [ij] could show that two differently 
transmitting channels contribute to transport. 

Theoretical studies of this phenomenon arc, on the one 
hand, based on calculations using density functional the- 
ory (DFT), [H, [H, [13, [S [H. In an early publication 



Wang and Berggren showed how Coulomb interaction can 
split the energy levels of up- and down electrons in a 
quasi one-dimensional system [ll|. They used DFT cal- 
culations with Hartree- and exchange potentials in local 
density approximation. Their findings were confirmed by 
more sophisticated calculations which include exchange- 
correlation potentials and take into account realistic gate 
potentials 17, 18, 19]. The observed difference of the up- 
and down energy levels gives rise to spin dependent trans- 
missions which manifest themselves in a 0.7 feature in the 
total conductance. However, to our knowledge there are 
no DFT results showing the correct temperature depen- 
dence. 

Besides DFT calculations, there are various theoreti- 
cal models describing different aspects of the 0.7 anomaly. 
Some models are based on the presence of spin-splitting 
assuming a density dependent separation 
of the up- and down energy levels from the beginning. 
These models can qualitatively reproduce the correct 
magnetic field and temperature dependence of the 0.7 
structure and are also suitable to describe shot noise [ij] . 
In a complementary approach the 0.7 anomaly is related 
to the Kondo effect [23l . [2^ by treating the quantum 
point contact as an interacting two-level system for the 
different spins. Qualitatively, this approach also leads to 
the observed temperature and magnetic field behavior of 
the 0.7 feature ^]. Very recently, also shot noise was cal- 
culated within this model showing agreement with exper- 
imental data [2^. Furthermore, interaction with phonons 
is used to explain the unusual temperature dependence 

In this work we present a comparatively simplified 
approach to the problem, which is based on the non- 
equilibrium Green function technique where the inter- 
action is incorporated at the Hartree-Fock level. We 
shall, furthermore, approximate the screened Coulomb 
interaction between the electrons by a repulsive contact 
potential. The fact that we can, within this approach, 
reproduce all relevant features of the 0.7 anomaly at 
temperatures close to zero, including the recently ob- 
served modification of the shot noise factor 14^, sup- 



2 



ports arguments in favor of tfie spin-splitting mechanism 
[11 [S 111 [H, [SO, HH, [2! (in line with experimental evi- 
dence provided in Ref. [13]) and indicates that the effect 
is rather robust with respect to the precise theoretical de- 
scription of the process. The short-range potential is fur- 
thermore chosen with regard to possible future transport 
experiments of ultracold fermionic atoms which precisely 
interact via the contact potential that we are using. 

This paper is organized as follows: in Section |TT| we in- 
troduce our model and present the relevant expressions 
that are used to calculate the transport properties. We 
show our numerical results in Section IlIII where we con- 
centrate on the influence of the coupling constant and 
the magnetic field on the conductance. We discuss the 
zero field case and show results concerning shot noise and 
finite temperatures. In Section HVl we summarize our re- 
sults and discuss transport of fermionic atoms through 
a constriction. The Appendix consists of a part about 
determining the strength of the interaction constant and 
other model parameters in an ideal two-dimensional elec- 
tron gas. Another part contains a detailed description of 
how to extend the recursive Green function algorithm to 
non-equilibrium processes. 

II. THE MODEL 

We describe a two-dimensional electron system with an 
additional in-plane confinement potential Vconf (a;, y) that 
defines the geometry of the quantum point contact. The 
in-plane magnetic field B = {B, 0, 0) oriented towards the 
transport direction gives rise to a Zeeman term only. For 
moderate magnetic fields the orbital contribution van- 
ishes with the choice A = (0, — i?z,0) for the vector po- 
tential and z = for the location of the two-dimensional 
electron gas (2DEG). Therefore, the non-interacting part 
of the Hamiltonian of the system can be written as 

= ^^:^ + V,ontix,y) + 9iiBBa, (1) 

where m is the effective mass and a — ±1/2 is the spin 
quantum number. The spin-up and spin-down energy 
levels are separated by the Zeeman energy Ez = g^isB, 
where g is the effective gyro-magnetic ratio and the 
Bohr magneton. Within our model the interaction of two 
particles located at r and r ' is described by 

Vintir,r') = -/S{f~f'), (2) 

with interaction strength 7. This choice of the interac- 
tion can be interpreted as a simple model for an effi- 
ciently screened Coulomb potential. For a homogeneous 
2DEG the Thomas-Fermi screening length is of the or- 
der As ~ 5 nm and the width of a typical constriction 
is roughly W 20 nm, see Appendix |X| Therefore, 
we do not expect that our model gives an accurate de- 
scription of the interaction. But it provides a transpar- 
ent physical picture of the mechanism causing spin split- 
ting. Moreover, our Hamiltonian is particularly devised 



to predict transport features of neutral fermionic atoms, 
as discussed in section |IV| In that case it is a very good 
approximation to use delta-like interactions. 

The coupling constant 7 can be estimated by calculat- 
ing the total interaction energy for a screened Coulomb 
potential in Thomas-Fermi approximation, as done in 
Appendix [K\ We find that 7 ~ 27r x It? /{2m) gives a 
realistic order of magnitude for the interaction strength. 

To calculate the transport properties of the system 
we use the Keldysh Green function approach [28|, [29j . 
This approach is very general as it allows to treat in- 
teractions and to include finite temperatures and bias 
voltages. The physical properties are obtained from the 
retarded and lesser Green function, G"" and G^. The 
former can be used to calculate properties such as the 
conductance, see Eq. ([5]); from the latter we get the par- 
ticle density, see Eq. Within the Green function 
approach the interaction is treated in a self-consistent 
way and can be included via a proper self-energy. For 
our calculations we take into account the first order of 
the perturbation expansion. The corresponding retarded 
self-energy is usually written as a sum of the Hartree and 
Fock self-energies, and S^, and has the form [2^ 

S^(r,r') ^S{r-r')Y. j df" j dE V,^,{f,r") x 

xGf{f'\r'\E), (3) 
= ^ I dEV,,Ar,r')G^{r,r',E). 

In general, the Hartree self-energy is local, and it involves 
a sum over all spin directions, whereas the Fock self- 
energy is non-local and depends only on the lesser Green 
function of the same spin orientation. However, in our 
case of delta-interactions, Eq. both the Hartree and 
the Fock contribution are local. exactly compensates 
the (cr — (T)-tcrm of the spin sum in E^, and we easily 
obtain for the total interaction self-energy 

S-t(r) = S^ + S?,=7n_.(f). (4) 

Here, 

n„{?) = -^ j dEG<{r,r,E) (5) 

is the density of electrons with spin cr. 

The Hamilton operator of the interacting system, 
TC^ — -fE?^jj(r), is a sum of the non- interacting Hamil- 
tonian ID) and the interaction self-energy ([l]), which acts 
like an additional local potential. This potential is differ- 
ent for the different spin directions: a spin-up electron en- 
counters a potential which is proportional to the density 
of spin-down electrons, and vice versa. Hence, there is a 
repulsive interaction only between particles with opposite 
spin directions. Therefore, any small imbalance between 
the density of up- and down electrons is increased by this 
kind of interaction. 
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FIG. 1: The geometry of the system. Each point corresponds 
to one site of the numerical grid; the grey points indicate the 
region where the interaction is gradually switched on and off. 
The leads are coupled to the left and to the right as indicated 
by the white bars. 



To solve the transport problem we discretize the spa- 
tial coordinates. The derivatives in the Hamiltonian are 
then written as finite differences, and the Hamilton op- 
erator is represented by a block-diagonal matrix [30|. 
The diagonal matrix elements contain an on-site energy 
and all local potentials Hfj — 4?i^/(2ma^) -f V^onf + 
'^inti'^i) + gfisBa (with a the lattice constant). The off- 
diagonal matrix elements for neighboring sites i and j 
are Tiij = —h^/{2ma^) and zero otherwise. As shown in 
Fig. [l] the geometry of our system has the shape of a lin- 
ear constriction with a hard-wall confinement potential 
Ki:onf (a;, y) which is zero inside the scattering region, and 
infinite outside. 

Moreover, the system is coupled to semi-infinite leads 
that have the same width as the outer slices of the con- 
striction. The leads are in thermal equilibrium character- 
ized by the chemical potential /i, and there is no effective 
electron-electron interaction in the leads. The interac- 
tion potential is gradually switched on/off in the narrow- 
ing region indicated by the grey points in Fig. [TJ The 
coupling to the leads can be exactly taken into account 
by self-energies El and for the left and right lead, re- 
spectively m, [sil • With these ingredients it is possible 
to calculate the full retarded Green function by matrix 
inversion 

G;(i?) = [i?-i/o--Sr„t-SL-S«]-\ (6) 

where Hq is given in Eq. (jT]). The Green function is a 
matrix of dimension 2N x 2N, where N is the number 
of lattice sites. Hence it would be very time consum- 
ing to invert the complete matrix in one step, as the 
computing time scales like N^. However, it is possible 
to implement a recursive algorithm that calculates the 
Green function of single slices of the scattering region 
and couples the slices via a Dyson equation. The details 
of this algorithm are explained in Appendix |B1 The re- 
cursive scheme scales with the third power of the width 
and only linearly with the length of the system. Thereby 
it is much more efficient than a direct matrix inversion. 
From the retarded Green function we get the lesser 



function using the kinetic equation 

g<{e) = g:{e)j:<{e)gue), (7) 

where the advanced Green function is obtained by her- 
mitian conjugation, G° = [GJ^]^. The lesser self-energy is 
T.<{E) = -2i f{E, fi)lm{j:L + ^r), where f{E,fi) is the 
Fermi-Dirac function. This relation holds as the leads 
are assumed to be in thermal equilibrium. So the lesser 
self-energy can be interpreted as the in-scattering rate 
for particles with energy i5 at a chemical potential /i. 
The lesser Green function (E) determines the particle 
density and hence the interaction self-energy, according 
to Eqs. ([3]) and ([5]). Thus, the interaction self-energy can 
be calculated from the retarded Green function, but in 
turn the retarded Green function depends on the inter- 
action self-energy. Hence, Eqs. (|3]) and ^ have to be 
solved simultaneously. 

The solution is carried out in an iterative way: we start 
with an initial guess for the interaction self-energy to cal- 
culate the retarded Green function with Eq. From 
this we get the lesser Green function, Eq. ([7]), and com- 
bining Eqs. ^ and ([5]) we obtain a new value for the in- 
teraction self-energy. We continue with this scheme until 
we have reached convergence. As soon as we have found a 
self-consistent solution we can calculate the conductance 
of the system using the Landauer formula 

G=^Tr{rLG''r;jG'^}, (8) 
with Tl/r = i (Si/fl, - ^l/r) ■ 

III. NUMERICAL RESULTS 
A. Dependence on the coupling strength 

We first calculate the conductance of the previously 
described model for zero temperature. It is convenient 
to use as the energy unit Ei — n"^ / {2mW^) , the en- 
ergy of the first transverse mode in the narrow region 
of the scatterer of width W. To break the symmetry 
between electrons with different spins we apply a small 
magnetic field so that the Zeeman energy has a value 
Ez — giJ^sB = 0.0015 El. The case of zero magnetic 
field is discussed separately in section HH CI 

The conductance for different interaction strengths 7 
is shown in Fig. [21 We find that for a small coupling 
constant 7 = 3.7 x %^ /{2ra) the up- and down contri- 
bution G| and G| differ from each other. This differ- 
ence is not due to the Zeeman shift, as the Zeeman en- 
ergy is approximately two orders of magnitude lower. It 
is caused by the effective repulsive interaction between 
electrons with different spin orientations. If the interac- 
tion strength is increased the up- and down contributions 
split more and more. Additionally, a small shoulder de- 
velops in the curve for the total conductance at values be- 
tween 0.6 and 0.7 of the conductance quantum. This is in 
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FIG. 2: (Color online) Total conductance (solid line) and the 
up- (dashed) and down (dotted) contributions for different 
interaction constants 7. The coupling constant takes values 
7 = 3.7 (black), 7 = 4.1 (red), and 7 = 4.5 (blue) from left to 
right in units of Ti^ /{2m). The curves for 7 > 4.1 have been 
horizontally offset for clarity. 



agreement with experimental results for the 0.7 anomaly 
in Ref. 0]. For sufSciently high interaction constants 
the contribution of one spin component to the conduc- 
tance drops down while increasing the chemical potential. 
These spin-resolved conductance curves coincide with re- 
sults obtained from transmission across a saddle potential 
in the presence of a Gaussian spin-splitting Q , and with 
corresponding DFT results [17]. 

The curves in Fig. [2] already indicate that the spin- 
splitting has to be sufficiently strong in order to get a 
visible effect in the total conductance. This is in accor- 
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FIG. 3: Lower panel: conductance curves for different values 
of the coupling constant (7 is given in units of /{2m)). 
The derivative (in arbitrary units) of each curve is plotted 
straight above the particular graph. The dashed lines mark 
the conductance value of the "plateau" for the lowest value of 
7, where a local minimum in dG/dn is visible. The different 
curves are offset horizontally. 



dance with the experimental results in Ref. [12] where the 
authors measure spin resolved contributions to the total 
conductance of a point contact. They find that even in 
samples that are not exhibiting a 0.7 feature the spin-up 
and -down electrons contribute differently to the total 
conductance. Fig. [3] shows how the 0.7 plateau devel- 
ops upon increasing the coupling strength 7. The lower 
panel shows the conductance curves and the upper panel 
the corresponding derivatives dG/dfi. The derivatives 
change from a single peak to a double peak shape as the 
0.7 plateau develops. The second peak in the derivative 
appears at a coupling constant 7 « 3.9 x Ti^ /{2m); the 
corresponding value for the plateau is G « 0.76 Go- In- 
creasing the interaction constant, the plateau gets more 
and more pronounced and eventually converges towards 
Go/2. Hence, in our model the interaction parameter 7 
governs the position and the width of the 0.7 plateau. 
The higher the conductance value at the plateau, the 
smaller is its width. There is no plateau above 0.76 Go 
in our model. 

To get information about the energy levels of the dif- 
ferent spin orientations we calculate = (S^Jjt)- 
The energy levels of non-interacting electrons E^ = 
El -\- q^bBct are shifted by the average interaction po- 
tential (Sj^t) felt by a particle with spin a. The brack- 
ets (...) denote a spatial average over the region with 
full electron-electron interaction (black points in Fig.[T]). 
The resulting curves in Fig. 2] show that the energy lev- 
els E„ are located around Ei for /i = 0. They are very 
weakly split by the Zeeman energy. With increasing /i 
the levels rise in energy as the constriction is populated 
with electrons and then start to split distinctly, as soon 
as the chemical potential is comparable with the energy 
levels E„ . The reason is that due to the Zeeman splitting 
the down-level is populated already at a lower chemical 
potential causing an imbalance between the density of 
spin-up and spin-down carriers in the constriction. The 
repulsive interaction between opposite spins tends to in- 
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FIG. 4; (Color online) Energy levels _E| (black) and Ey (red) 
for 7 = 4.5 X h'^/{2m) and Ez = 0.0015 Ei. The dashed curve 
shows the chemical potential ji. 
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crease any imbalance. A small excess of down electrons 
repels up electrons from the constriction, which results 
in a larger excess of down particles. 

The spin-splitting vanishes when the chemical poten- 
tial is well above both energy levels. In the range of the 
chemical potential where the up- and down energy levels 
are split, also the up- and down contribution to the con- 
ductance differs, as shown in Fig. O The obtained energy 
levels shown in Fig. |4] are in line with DFT results IITI. 

Before comparing with the spin-splitting models |2d . 
[2H we shall note that the quantities plotted in Fig. [3] 
are only estimations for the energy levels. Due to the 
geometry of our system the transverse modes are broad- 
ened with a width of the order of i?i/2, as can be seen 
in the conductance curves of Fig. [2l Therefore, our re- 
sults seem to confirm the assumption of the spin-splitting 
models, that the energy levels start to split as soon as the 
chemical potential crosses the up- and down energy lev- 
els. Additionally, we observe a "pinning" of the upper 
energy level to the chemical potential within a substan- 
tial range of /x, that means E-^ evolves parallel to the 
chemical potential /i right after the splitting. The pres- 
ence of this level pinning is essential in the spin-splitting 
models in order to get a 0.7-plateau. In our calculations 
the plateau also appears in the range where level pinning 
is present. 




chemical potential [i I 

FIG. 5: Conductance curves for 7 = 4.5 x 'h? /{2m) and differ- 
ent magnetic fields. The corresponding Zeeman energies vary 
from to 0.14 Ei in steps of 0.012 Ei from left to riglit. The 
curves are horizontally offset. 



value is lower than in experiments where fields of about 
10 T are necessary to get a plateau at 0.5 Go 0]. 



C. The zero-field case 



B. Magnetic field dependence 

The shape of the conductance curves is influenced by 
the magnetic field. Fig. O shows that the 0.7 plateau 
evolves from a small shoulder at G ~ 0.65 Go to a wide 
plateau at G = 0.5 Go as the magnetic field is increased. 
This is in agreement with experiments [T^ . In the 
high- field limit Zeeman splitting is the dominant effect. 
The energy levels of the different spins are separated by 
the Zeeman energy which causes a plateau at 0.5 Go 
even in the case of non-interacting electrons. The rea- 
son is that spin-down electrons contribute to transport 
at chemical potentials fi ^ Ei — g^sB, whereas for spin- 
up electrons /x > Ei+gfj,BB has to be fulfilled. For strong 
magnetic fields the effect of electron-electron interaction 
is only to broaden the Zeeman spin-split plateau at one 
half of the conductance quantum. 

For a more quantitative comparison between our re- 
sults and experimental data it is useful to re-scale our 
quantities and give the magnetic fields in units of Tesla. 
Therefore, we have to associate an energy value with 
El = h^TT^ /{2mW^). If we insert the approximate width 
W = •\/2/(7rn) of a quantum point contact from Eq. (|A2[) 
we get 

„ 1 h-K^n Ez , , 

B = Ez ~ -. (9) 

gfiB 2ge Ei 

The maximum field applied in Fig.[5]then corresponds to 
B « 5.8 T where we used g = 0.44 for bulk GaAs and a 
density of n = 1.8 x 10^^ cm~^ [5]. This magnetic field 



For the previous calculations we always applied a fi- 
nite magnetic field. Due to this field the energy levels of 
the electrons with different spins were separated so that 
the down state can be populated at smaller chemical po- 
tentials than the up state. That is the reason why up 
electrons are repelled from entering the constriction, as 
down electrons are already present at a lower chemical 
potential. So the repulsive interaction between particles 
with opposite spins leads to an enhancement of an ini- 
tially small asymmetry between the density of up- and 
down electrons in the constriction. 

However, in the case of zero magnetic field the Hamil- 
tonian ^ together with the interaction self-energy (|4]) 
is strictly symmetric with respect to spin-up and spin- 
down. Therefore, the resulting conductance curves also 
have to show the same symmetry. The result for B — 
and 7 = 4.5 x fi^ / {2m) is displayed as the dashed curves 
in Fig. [HI As expected the contributions of the up- and 
down electrons exactly coincide and the total conduc- 
tance has no additional features below the first step. 

We can investigate the stability of the symmetric solu- 
tion by slightly disturbing the symmetry of the system. 
For each point calculated we start with a small magnetic 
field E2; — 0.0015 El. During the first four steps of the 
self consistency loop we turn off the magnetic field ac- 
cording to E'^ / s where s is the number of the iteration 
step. After four steps we set the magnetic field exactly 
to zero and continue iterating until the results are con- 
verged. In that way we obtain an asymmetric solution 
for the spin-up and spin-down contributions. In contrast 
to the finite field results depicted in Fig. [21 here the split- 
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FIG. 6: (Color online) Conductance in the zero-field case. 
The dashed lines are for B — 0; there is no difference between 
the t and ^-contribution and the total conductance does not 
exhibit a 0.7 feature. The solid lines show the results for 
B where the t (red) and j (green) contributions differ 
and the total conductance (bold blue line) has a shoulder at 
G fa 0.65 Go. 



ting sets in abruptly at a chemical potential ii — 1.20 i^i. 
In the range where G| and Gi are different a shoulder 
appears in the total conductance G. Those points where 
the spin-splitting is absent coincide with the points for 
B — 0. So the symmetric solution with = G| is un- 
stable and we find a 0.7 anomaly even in the case of zero 
magnetic field. 

In our case the down-contribution to the conductance 
dominates when we apply a positive magnetic field. With 
a negative field the different spin directions would change 
their roles. In reality the asymmetry between spin-up 
and spin-down may be caused by residual magnetic fields 
or temporal current fluctuations. Also magnetic impuri- 
ties, as well as nuclear spins and dynamic nuclear po- 
larization might play a role in breaking the up- and 
down-symmetry. Our numerical results show that a very 
weak asymmetry is sufflcient to get spin-splitting. We 
obtained spin-split results for Zeeman energies down to 
Ez = 3 X 10^^ El, corresponding to a magnetic field 



strength of about B w lO"'^ x. 



D. Shot noise 

In recent experiments shot noise was measured in 
quantum point contacts exhibiting a 0.7 anomaly [T3l.[l^. 
In the framework of Landauer-Biittiker theory the shot 
noise power S" in a two-terminal device is given by [33 | 



h ^ 



dET„, 



Here, fL/R is the Fermi distribution function of the 
left /right contact. If the energy scale on which the trans- 



mission functions Tn,a{E) vary is large compared to tem- 
perature kgT and applied source-drain voltage eVsd, the 
transmissions can be treated as constants. Then the en- 
ergy integral over the distribution functions can be per- 
formed, yielding 



S = 2J\f^ 



2e2 



eVsd coth 



- 2kBT 



with the noise factor defined as 



(10) 



(11) 



The noise factor of one single channel vanishes for zero or 
perfect transmission, and it is maximal for Tn^a = 1/2. 

By simultaneous noise and conductance measurements 
it is possible to extract information about spin-resolved 
transmission coefficients T| and . Whereas the conduc- 
tance is proportional to the total transmission, Ttot = 
+ Tl, the noise factor in the single mode case is 
TV = 5Ttot(l - Ttot) + T^T^. Only in the case of non- 
interacting particles where T| = T| = T the noise factor 
reduces to TV = T(l -TV 

The authors of Ref. lj| measured the shot noise power 
and fitted their experimental results with Eq. ()10|1 using 
TV as fitting parameter. They find a suppression of noise 
around the anomalous conductance plateau. That gives 
experimental evidence that near the 0.7 feature electrons 
are transported by two channels with different transmis- 
sions, as also stated in [l^. This agrees with our re- 
sults for G| and G^ displayed in Fig. [21 For conduc- 
tance values between and 1 the experimentalists find 
an asymmetric dome shape for the noise factor evolv- 
ing into a symmetric double-dome structure by applying 
a magnetic field. They are also able to reproduce this 
behavior with Reilly's phenomenological model [2l[. In 
a recent publication the same noise characteristics was 
obtained using a Kondo model 
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FIG. 7: (Color online) Noise factor for 7 = 4.0 x h? /(2m) and 
for different values of the magnetic field Ez/ Ei = 0.0, 0.0007, 
0.006, 0.012, 0.018, 0.23 from top to bottom. The dashed line 
indicates the result for non-interacting electrons. 
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The noise factor within our model is depicted in Fig.[7| 
for 7 = 4.0 xh^ / {2m) and for different magnetic fields. In 
agreement with the results in Ref. [l^ we find an asym- 
metry of the noise factor with respect to G = 0.5 Go. The 
shot noise is suppressed at conductance values around 
0.7 Go which accounts for the differently transmitting 
channels in that range. The spin-down channel has al- 
most perfect transparency and hence does not contribute 
to the noise factor whereas for G ~ 0.3 Go both channels 
are equally transmitting and contribute equally to the 
noise factor. 

By increasing the magnetic field a second maximum 
appears and the noise factor evolves towards asymmet- 
ric shape. In contrast to the model results in [ij] where 
the maximum of the right dome is stationary it first 
rises slightly in our model and then drops down again 
for Ez/Ei > 0.04. However, for very strong magnetic 
fields, Ez/Ei = 0.23, {B « 9.5 T), the noise factor is 
symmetric with two maxima at Af = 1/8. This accounts 
for spin resolved transmission of electrons due to the Zee- 
man splitting. The discontinuity of the lowest curve at 
G « Go/2 is caused by the small oscillations of the con- 
ductance around G — Gq, see e.g. Fig. [21 In that regime 
one finds two different noise values for one conductance 
value. 



E. Temperature dependence 

The 0.7 anomaly is accompanied by a peculiar temper- 
ature dependence: within a certain range the 0.7 plateau 
gets more pronounced if the temperature is increased 
This behavior can be reproduced by the spin- 
splitting models [l^, |2l[ as well as by the Kondo model 
[23| and by interaction with phonons 25]. However, to 
our knowledge there are no DFT results exhibiting such 
a temperature behavior p7| . 

If we include finite temperatures in our calculations we 
find a reduction of the spin-splitting as shown in Fig. [S] 
The difference between the transmission of up and down 
electrons is reduced with increasing temperature which 
makes the 0.7 plateau less pronounced. This result con- 
tradicts the experimental findings. As DFT calculations 
are also not able to capture this phenomenon it is possible 
that a mean-field description is not sufficient to explain 
the temperature dependence. 

In contrast to our approach the spin-splitting models 
qualitatively yield the correct temperature dependence 
[20l . [2l|. In the model approach, however, the temper- 
ature affects only the computation of the conductance 
whereas the spin-splitting is assumed to be temperature 
independent. In our calculation the temperature also en- 
ters in the computation of the density, Eq. ([5]), as G< (E) 
is truncated around the Fermi level plus several fc^T. 
Hence, the densities and thus the interaction potentials 
depend on the temperature which results in a temper- 
ature dependent spin-splitting. The spin-gap vanishes 
at temperatures ksT > 0.08 Ei. Second, the spin- 




chemical potential ji / 



FIG. 8: (Color online) The conductance contributions G| 
and Gi for 7 = 4.5 x /{2m) and for different temperatures 
ksT = (solid black line), ksT = 0.029 Ei (dashed red line), 
and ksT — 0.058 i?i (dotted blue line). The inset shows the 
corresponding total conductance G = + G^. 



splitting models assume sharp energy levels with step- 
like transmission functions Ta{E) — Q{E ~ Ea-), where 
Q{x) is the Heaviside step-function. Finite temperatures 
lead to a smearing of the conductance and a 0.7 struc- 
ture is found for teinperatures smaller than the spin-gap, 
ksT < \E-^ — Ei\ [20]. In our model the energy levels 
exhibit a broadening due to the geometry of the system 
even at zero temperature, ksT — 0. The broadening is of 
the order Ei/2 (see Fig. [5]), larger than the level splitting 
(see Fig. [3|). Hence, allowing for finite temperatures the 
broadening is further enhanced which leads to a decrease 
of the 0.7 plateau. 



IV. DISCUSSION AND OUTLOOK 
A. Electron transport 

The presented model describes transport of locally in- 
teracting electrons. In Hartree-Fock approximation only 
electrons with different spins are interacting repulsively. 
This suggests a very intuitive physical picture: if the 
scattering region is predominantly occupied by one spin 
species, electrons with opposite spin are repelled from the 
constriction. So this kind of interaction favors an asym- 
metric population of the quantum point contact. Despite 
its simplicity the model is adequate to qualitatively ex- 
plain different aspects of the 0.7 anomaly. We see how 
interaction can cause an asymmetry between the spin-up 
and spin-down transmission resulting in a shoulder in the 
total transmission. The magnetic field dependence of the 
0.7 feature is well reproduced and we find an instability 
phenomenon in the zero field case leading to spontaneous 
spin polarization. Our model also accounts for shot noise 
suppression at the 0.7 plateau. 

Going beyond Hartree-Fock one expects that spin- 
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splitting is weakened or even vanishing in the strict one- 
dimensional case. However, it was shown by exact meth- 
ods that Hubbard chain models can have a ferromag- 
netic ground state if one is not restricted to exactly 
one-dimensional systems Is^l , in accordance with the 
Lieb-Mattis theorem [sHi- So in our system which is 
based on a two-dimensional description a spin polarized 
ground state is possible. 

Because the assumption of a delta-like interaction po- 
tential might be too crude for the junction in the 2DEG, 
we also performed calculations with Coulomb interaction 
and full exchange. For the conductance curves we found 
similar results as for delta interaction, as shown in Fig. [9] 
for different coupling strengths. The dimensionless cou- 
pling constant can be estimated as 



7c 



ATTEoe h 



0.20, 



(12) 



where m = 0.07 mg, e = 13 and a — Inm was inserted. 
Replacing the Coulomb interaction by a Yukawa poten- 
tial V{r) ~ exp(— r/A)/?" we find that the results evolve 
towards the curves for delta interaction with decreasing 
screening length A. So the effect of spin-splitting remains 
robust even in the limit without screening. In the case of 
Coulomb interaction the diagonal part of the Fock self- 
energy also compensates the short-range contribution to 
the Hartree self-energy, Eq. ([3]). This leads to an effec- 
tive short-range repulsion between different spins, sim- 
ilar to the case of delta-interaction. We therefore can 
conclude that the repulsive interaction between electrons 
with opposite spin causes spin-splitting even in the case 
of long-range Coulomb interaction. 

Our approach cannot reproduce the experimentally 
observed temperature dependence of the plateau struc- 




2 3 4 

chemical potential n / 

FIG. 9: (Color online) Total conductance (solid line) and 
the up- (dashed) and down (dotted) contributions for differ- 
ent interaction constants 7c calculated with full Coulomb in- 
teraction. The dimensionless coupling constant takes values 
7c = 0.18 (black), jc = 0.20 (red), and 7c = 0.22 (blue) from 
left to right. The curves for 7c > 0.20 have been horizontally 
offset for clarity. 



ture |5|, which is also the case for DFT calculations 
[TgI UJ\ . [TsI . [T9| . This admits two possible interpreta- 
tions: On the one hand, Kondo-type correlations could 
be responsible for the temperature-induced enhancement 
of the 0.7 feature. This mechanism was theoretically 
suggested in Ref. (23j and experimentally supported in 
Ref. |10j , but still awaits ultimate confirmation by ab ini- 
tio calculations that are able to take into account such 
correlations and do not involve any tunable parameter. 
It was, on the other hand, suggested [l^ that phonons 
could be at the origin of this effect. 



B. Transport of fermionic atoms 

To discriminate between these two complementary in- 
terpretations, we propose to perform transport experi- 
ments with uUracold fermionic atoms, such as ®Li for in- 
stance, which can nowadays be routinely confined within 
magnetic or optical trapping potentials and cooled down 
to temperatures close to the BCS transition ^3^. In the 
context of interaction-induced modifications of the con- 
ductance, optical (rather than magnetic) techniques for 
the confinement of the atoms would be required in or- 
der to trap both spin species of the fermionic atom. A 
quasi two-dimensional configuration could, for instance, 
be realized by a rather strong one-dimensional optical 
lattice which creates a sequence of disk-like confinement 
geometries for the atoms, and a matter- wave guide with a 
constriction could be induced by additional laser beams 
that are focused onto the disk within which the atoms 
are confined. 

According to Ref. [sj, the effective interaction con- 
stant 7 that characterizes the contact potential ^ would, 
in the case of two-dimensional ultracold fermions, be 
given by 



1 



27r^ -f ln(^ 

as V ttE 



(13) 



Here, m is the mass of the atom, luj_ denotes the fre- 
quency of the harmonic confinement in the transverse 
direction (i.e., along the "third" dimension), aj_ = 
yJhfnvjjZ is the corresponding oscillator length, E de- 
notes the total energy of the collision process between 
two atoms in the center-of-mass frame, and represents 
the s-wave scattering length between two atoms with op- 
posite spin. Both length scales, Og and aj^, can be manip- 
ulated, via Feshbach tuning (see, e.g., Ref. 38]) as well 
as through the intensity of the optical lattice. It would 
therefore be possible to realize configurations for which 
the effective interaction strength 7 is of the order of the 
values that were discussed in Section UlIl 

To measure the atomic 0.7 anomaly, we propose to pre- 
pare the fermionic atoms in a large double- well trap that 
is optically created within the two-dimensional confine- 
ment geometry, and let them escape from one well to the 
other through a small "bottleneck" corresponding to the 
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constriction of Fig. [TJ Counting the number of atoms 
that are transported across the bottleneck within a finite 
time scale should give rise to a current of atoms close 
to the Fermi level. This current can be directly trans- 
lated into an "atomic conductance" in a similar way as 
in Ref. [3§| , which would also display a step-like behavior 
when the height of the constriction is lowered by optical 
techniques. Magnetic fields can again be used to break 
the symmetry between spin-up and spin-down fermions, 
and the temperature could possibly be controlled by let- 
ting the fermionic cloud interact with a gas or condensate 
of bosonic atoms (e.g., by preparing a mixture of ^Li and 
^Li atoms). As phonons are clearly absent in this setup, 
any observed feature in the 0.7 anomaly that is not re- 
producible by mean-field approaches would necessarily 
be due to (Kondo-type) correlations. 

In short summary, is should be possible to realize 
transport experiments with ultracold fermionic atoms 
where the 0.7 anomaly in the conductance would be ob- 
served. We expect that such experiments would provide 
new insight into the central mechanism that underlies 
this phenomenon. 



When the first channel opens the effective width can be 
estimated by equating the chemical potential with the 
energy of the first sub-band fiLLj/2 for a parabolic con- 
finement. The width of the confinement potential at this 
energy is = l^hj (rnui), which gives 



(A2) 



W = ^ — 



So we find that the effective width of a quantum point 
contact is of the order W — 20 nm. Inside the constric- 
tion the density is expected to be lower than in the ho- 
mogeneous 2DEG, so the effective width will be larger 
than the above estimated value. 

For (5-interaction the coupling constant is given by the 
spatial integral over the interaction Hamiltonian. So 
we calculate the corresponding quantity for a screened 
Coulomb potential 



d r - 



2eQe' 



(A3) 



Inserting the screening length, Eq. (jAip . we find 
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APPENDIX A: ESTIMATION OF THE MODEL 
PARAMETERS 

The coupling constant 7 is the main parameter of our 
model. It has to be sufficiently high in order to get an ob- 
servable effect of the electron-electron interaction. Here 
we want to estimate an upper limit of the interaction 
strength using an exponentially screened Coulomb po- 
tential. 

In a homogeneous 2DEG the screenin g le ngth in 
Thomas-Fermi approximation is given by [40| 



7 



27r- 



2m 



(A4) 



A. 



2-Keoefi 



(Al) 



where e is the average dielectric constant of the two ma- 
terials on both sides of the 2DEG. For a GaAs/AlGaAs 
interface we find As ~ 5 nm, where m — 0.07 mo and 
e = 13 was used. To compare the screening length with 
the width of the constriction, we have to estimate the 
typical dimensions of a point contact. The lithographic 
width is of the order of several hundred nanometers, but 
the electrons are confined by the electrostatic potential 
due to the gates. The effective width of the constriction is 
then controlled by the gate voltage. From experiments we 
know the typical density of carriers n = 1.8 x 10"'^^ cm~^ 
which is related to the chemical potential /i = h^Tm/m. 



This is just a rough estimation as several aspects are 
neglected. First, in Thomas-Fermi approximation the 
screening length in two dimensions is independent of the 
electron density. But beyond this approximation one 
finds an increasing screening length as the charge den- 
sity goes to zero [40]. This reflects that screening is less 
efficient if the particle density is too small. 

Second, screening in two dimensions is not as strong as 
in three-dimensional systems. The asymptotic behavior 
of the screened potential is not exponential, but it follows 
an r~^ law [40| . However, the resulting coupling constant 
does not differ dramatically from the one obtained by 
exponential screening. Both facts would give rise to an 
even higher upper limit of the coupling strength. 



APPENDIX B: RECURSIVE ALGORITHM FOR 
NON-EQUILIBRIUM GREEN FUNCTIONS 

The recursive Green function algorithm is widely used 
for calculating electronic properties of two- and three- 
dimensional systems. The basic idea is to build up the 
full Green function slice by slice instead of evaluating it in 
one step. Thus, the dimensions of the matrices that have 
to be inverted are strongly reduced. If the Green function 
Go of a semi- infinite region and an adjacent isolated slice 
is known, it is possible to calculate the Green function G 
of the coupled system using Dyson's equation 



G = Go + GqVG. 



(Bl) 



(For this derivation we omit the spin index a and the su- 
perscript r for the retarded functions). Here, V denotes 
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n n+1 N+1 

FIG. 10: The Green function is constructed by coupling single 
slices starting from one of the leads (grey regions). 



the hopping matrix between the two adjacent shoes. The 
Green function of a semi-infinite lead can be calculated 
analytically [3l[. So it is possible to start with an iso- 
lated lead and then add slice by slice until the opposite 
lead is reached. This is schematically shown in Fig. [TOl 
After coupling that lead to the rest of the system one has 
obtained the Green function of the complete system at 
the surface of one lead. This Green function contains all 
information to calculate the current through the system. 
The above described procedure is explained for example 
in Refs. ^\^. 

As we are interested also in the electronic density which 
is determined by the lesser function G^, the above ex- 
plained algorithm is not sufficient. Here we present an 
extension to the usual recursive Green function method 
which allows us to calculate the retarded Green function 
between the two leads as well as the lesser Green function 
(see also Ref. [1^). The condition to apply this algorithm 
is that all relevant self-energies are diagonal so that the 
effective Hamiltonian that has to be inverted, Eq. (O, 
keeps its block-diagonal structure. This condition is ful- 
filled for the Hartree self-energy and also for the Fock 
self-energy in our case of delta interaction. In the gen- 
eral case of full Coulomb interaction the Fock self-energy 
is not diagonal, see Eq. ([3]), so the presented method can 
not be used. 

We first show how to add one single slice to a semi- 
infinite region. In the following we use the notation G'^'-"-' 
for the Green function of an isolated slice n, and G^/^*^"^ 
for the Green function of the right/left semi-infinite re- 
gion starting at slice n. The full Green function of the 
complete (infinite) system is denoted by G (without su- 
perscripts). In order to couple the Green function G'^'^"-' 
of the isolated slice n to the Green function G''^*^""''^^ that 
covers all lattice sites to the right of (n+1), we use the 
Dyson equation (|Bip . 



Qliin) 



n ) = { n 



^QSin) QB.(n+l) 

(n|(G^(")+G^("+i))l/G«(")|n). 



As G^'^"+^' has no matrix elements with slice n, the 
terms (n 

|G^("+i)|n> vanish and we get 



n ) = { n 



n) + 



^(n|G^("V>(a|^l&)(&|G^^" 



Noting that G'^^") has only non-zero matrix elements 
with slice n we get the constraint a — n. As the cou- 
pling matrix V acts only between adjacent slices and has 
no overlap with other slices b is restricted to the values 



,± 1. With 



1|G«(' 



n)| 



1) = we find 



(B2) 



where Gn.m — {n\G\m) is the sub-matrix of G related to 
the slices n and m. The Green function G^^"''„ appearing 
in Eq. (|B2p can be calculated via the Dyson equation in 
a similar way, and we get 



G 



R{n) _ QR{n+l 



n+1, n+1 ^ n+l,n^n,n 



^i?(n 



Inserting this result into Eq. (|B2[) and solving for G 
we obtain 



R(n) 



QR(n) _ 



[E - Hn,n) - K.n+lG, 



i?(n+l) 
n,n + l'~^n+l,n+ 



iK 



+ l,n 



where we used Gn^n^ 



(B3) 

= {E - Hn,ny^. Therefore, 
Eq. (|B3p allows us to calculate the Green function cover- 
ing all lattice sites to the right of slice n from the Green 
function to the right of (n+l). In that way we have added 
one slice. Iterating this scheme we can finally obtain the 
Green function G^ \ at the left end of the scattering re- 
gion. Then one has to connect the Green functions of the 
two semi-infinite sections to get the full Green function 
Gi^i (without superscript) at the left end of the scatterer. 
This we obtain by using the Dyson equation 



(1|G|1) 



gL(0)^gi^(l)- 

l|(G^(°) + G^«)yG|l), 



and we find 
Gi.i — 



In this equation Gq g"'* is the surface Green function of the 
semi-infinite left lead. The Green function Gi^i contains 
all information about the refiection coefficients at the left 
lead. 

In an analog way we can start from the left lead and 
calculate all Green functions from left to right by 



QL(n) 



{E - Hn,n) - ^n,n-l'-'n-"l,n-l^n 



-l,n 



(B5) 

and finally obtain the full Green function at the right end 
of the scatterer, 



G 



N,N 



1 r^Wv 



N,N ^N,N+l'^N+l,N+l^ N+1,N 



RiN+1) 



G 



L{N) 
N.N ■ 

(B6) 

Here G^^_^''^^_^j^ is the surface Green function of the right 



lead. 



11 



Knowing the full Green functions at both ends of the 
scatterer, Gi^i and Gjv.at, we can now compute the full 
Green functions between the ends and any slice n inside 
the scattering region. We use the Dyson equation 



{n\G\l) = 



'G 



L(n-l) 



G 



R{n] 



QL(n-l) ^ QR{n)\ yQ ^ 



to obtain 



Gn,l — G^^^'^Vn.n-lGn-l.l, 



(B7) 



where the Gn^'n^ 
gously one finds 



are calculated from Eq. (|B3p . Analo- 



GnM ~ G^}-'^VnM+\Gn+l,N 



(B8) 



with the Green functions G^J^ from Eq. (jBSp . The last 
two equations allow us to compute the Green functions 
G„^i and G„_Ar recursively by starting with the Green 
functions Gi^i and Gat^at at the ends of the scattering 
region. 

Now it is possible to compute the diagonal elements of 
the lesser Green function which are needed to calculate 
the electron density, Eq. ([5]). A diagonal matrix element 
of G"^ reads according to Eq. ([7]) 



[G" 



X2 



(B9) 



with i, j € {n = l,n = TV}. The self-energy is only 
non-zero at the ends of the scatterer where the lattice 
sites are coupled to the leads. So the indices i and j 
are from the first and last slice of the scattering region. 



Therefore, the Green functions calculated from Eqs. (|B7p 
and (|B8p enter here. 

The complete recursive procedure can be summarized 
in the following steps: 

• calculate and store all G^\^^ from the lead Green 
function G^^^^j^^ by means of Eq. (|B3p : 



• calculate and store all Gn^n' from the lead Green 

-,L(0) 



function Gq q by means of Eq. (|B5P ; 

• compute the full Green functions Gi.i and Gn.n 
at the left and right end of the scatterer using 
Eqs. dm and (|B6)) : 



• use Eqs. (IB7p and (jBSp to calculate and store all 
Gn,i and Gn,N from Gi^i and Gn,n- One of those 
Green functions is G]\r,i which contains information 
about the transport properties; 

• finally, one obtains the diagonal elements of G"^ 
from Gn,i and Gn,N by means of Eq. (jB9p . 

In total, one has to run four times through the entire 
system in order to be able to calculate G^ as well as parts 
of G^ which are needed for the reflection and transmis- 
sion coefficients. If the calculation of G^ is not necessary 
it is enough to pass the system twice to get all reflection 
and transmission coefficients. So the scheme reduces to 
the standard recursive algorithm [3l|, |4l| . If one is only 
interested in the current, passing the scatterer once is 
sufficient. After computing Gi^i with Eq. (|B4p the to- 
tal reflection is known. Employing current conservation 
(unitarity) it is possible to get the total transmission and 
hence the current. 
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